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1.  INTRODUCTION 


This  final  report  presents  the  results  of  a  12-month  study  of  adaptive 
filters  for  AMTI  radar.  The  results  are  applicable  to  several  types  of 
radars,  including  .il  and  AEW  systems.  Adaptive  systems,  and  adaptive  filters 
in  particular,  sense  the  existing  noise  field  and  optimize  a  set  of  system 
parameters.,  In  an  adaptive  filter,  the  filter  weights  and  transfer  functions 
are  adjusted  to  maximize  the  output  signal-to-clutter  (plus  receiver  noise) 
ratio. 

The  clutter  spectrum  in  an  AMTI  radar  is  a  function  of  several  variables 
Including  scan  angle,  radar  velocity,  antenna  pattern,  and  angular  distribu¬ 
tion  of  clutter.  Rain  backscatter  is  an  important  component  of  clutter  in 
higher  frequency  radars.  Its  location,  mean  radial  velocity,  and  velocity 
spread  due  to  wind  shear  are  not  generally  known  a  priori.  An  adaptive 
filter  senses  all  of  these  clutter  properties  and  optimizes  the  filter  re¬ 
sponse  for  detection  of  targets  with  a  selected  radial  velocity.  In  most 
cases,  a  bank  of  adaptive  filters  will  be  required  to  cover  the  target 
doppler  spectrum.  The  theory  of  adaptive  filtering  is  discussed  in  Section  2 
of  this  report. 

During  this  study,  a  computer  program  was  developed  for  the  investigation 
of  both  the  steady-state  response  and  transient  performance  of  adaptive  AMTI 
filters.  Provisions  for  both  ground  clutter  and  rain  clutter  were  included. 

A  variety  of  parameters  are  variables  in  this  program,  including:  radar 
velocity,  scan  angle,  target  radial  velocity,  adaptive  loop  parameters,  number 
of  pulses  processed  coherently,  and  rain  clutter  properties  (magnitude  relative 
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to  ground  clutter,  mean  radial  velocity,  and  spectral  width).  A  large 
number  of  cample  cases  were  run  to  develop  some  insight  into  the  performance 
of  adaptive  filters  and  to  uncover  problem  areas.  Selected  results  are  pre¬ 
sented  in  Section  3  of  this  report. 

The  suitability  of  adaptive  filters  for  future  AMTI  radars  depends  both 
on  the  performance  of  adaptive  systems  and  on  their  complexity.  Methods  of 
reducing  the  complexity  of  adaptive  control  loops  for  filters  (or  adaptive 
array  antennas)  were  also  investigated  during  this  study.  A  technique  for 
reducing  significantly  the  complexity  of  digital  adaptive  filters  is  discussed 
in  Section  4  and  in  Appendix  B  of  this  report. 

Suggested  areas  for  further  research  on  adaptive  filters  and  conclusions 
of  this  study  are  presented  in  Section  5.  A  simplified  and  more  accurate 
method  of  computing  the  control  loop  noise  power  in  adaptive  filters  is  con¬ 
tained  in  Appendix  C. 
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2.  ADAPTIVE  FILTERS 

Consider  an  airborne  coherent  pulsed  radar  designed  to  detect  moving 
targets  in  a  clutter  background.  Let  V  denote  a  column  vector  of  the  consecu¬ 
tive  received  signals  from  one  range  cell,  V  -  (v,,v,,,..,v„),  where  v  is 

x  X  4  tv  n 

the  complex  video  for  the  n^  sample  and  T  denotes  the  transpose.  An  11TI 
filter  output  is  obtained  by  forming  the  weighted  sum  of  these  v  ,  W^V,  where 
W  is  a  column  vector  of  complex  weights  w^.  The  filter  response  is  determined 
by  the  w^,  which  are  chosen  to  maximize  the  output  signal-to-interference 
ratio  (S/C). 

The  interference  (clutter  plus  receiver  noise)  power  in  the  output  is 

c  =  eJjwtv|2J  =  wJmw  (1) 

where  M  is  the  covariance  matrix  of  the  interference  process  with  elements 


Both  clutter  and  receiver  noise  (but  not  signal)  are  included  in  the  computa¬ 
tion  of  the  covariance  matrix.  In  these  equations,  *  denotes  the  complex 
conjugate  and  E  the  expectation  or  average.  The  output  voltage  samples  v^ 
and  weights  w^  are  complex  quantities  retaining  both  phase  and  amplitude 
information. 

In  designing  a  filter,  some  assumption  must  be  made  a  priori  concerning 
the  target  doppler  frequency.  Consider  the  case  where  a  separate  filter  is 
synthesized  for  each  of  a  set  of  target  doppler  frequencies,  analogous  to  a 
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filter  bank.  Let  S  denote  a  column  vector  of  signal  phases  for  the  K  con¬ 
secutive  samples,  =  (l.e^^e*^, . . .  ,e^K  A  signal  of  amplitude  A 

received  from  a  target  with  doppler  frequency  f  is  represented  by  AS,  with 
=  2T[fT,  where  T  is  the  pulse  repetition  period.  The  corresponding  signal 
voltage  at  the  filter  output  is  AW^S  and  the  output  S/C  ratio  is 


(3) 


In  most  cases  of  interest,  the  interference  processes  are  Gaussian. 

The  optimum  detection  algorithm  for  these  cases  is  to  form  a  maximum  S/C 
filter  and  compare  the  signal  amplitude  at  this  filter  output  with  a  detec¬ 
tion  threshold  The  filter  which  maximizes  S/C  has  weights  propor¬ 

tional  to 

W  =  M-1S*  (A) 


An  adaptive  filter  will  now  be  described  which  senses  the  interference  spectrum 
or  covariance  matrix  M  and  generates  a  weight  vector  that  approaches  Eq.  (A). 

A  K-pulse  adaptive  filter  is  implemented  with  K  separate  adaptive  control 
loops,  each  generating  one  complex  weight  wr,  as  illustrated  in  Figure  1.  To 
cover  all  range  cells  in  a  pulsed  radar,  K-l  delay  lines  are  required,  each 
with  a  delay  equal  to  the  pulse  repetition  period.  The  operations  indicated 
in  Figure  1,  and  in  the  control  loops  illustrated  in  Figure  2,  can  be  performed 

with  either  analog  or  digital  circuitry.  The  control  loops  are  identical 

A 

except  for  the  steering  signals  Sr,  which  are  matched  to  the  desired  signal 
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Figure  2.  Adaptive  Control  Loop 

doppler  frequency.  With  low-pass  filters  employed  as  integrators  in 
Figure  1,  satisfying  the  equation 

TU  +  u  «=  z  (5) 

n  n  n 

[2] 

the  weight  vector  satisfies  the  equation 

£  W  +  01  +  1/6) W  -  S*  (6) 

where  I  is  the  identity  matrix.  When  the  gain  G  is  large  (relative  to  the 
eigenvalues  of  M),  the  weights  approach  the  steady-state  solution  of  Eq.  (4). 

The  theory  of  adaptive  filters  is  closely  analogous  to  the  theory  of 
adaptive  array  antennas,  which  is  detailed  in  [2],  Since  the  covariance  matrix 
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M  is  Herraitian,  the  transient  response  of  the  system  can  be  computed  in 
normal  coordinates.  The  effect  of  control  loop  noise  in  these  systems  has 

r  2 1 

also  been  analyzed.  A  simplified  and  more  exact  derivation  of  the  con¬ 
trol  loop  noise  equation  is  contained  in  Appendix  C. 
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3.  AMT I  PERFORMANCE  OF  ADAPTIVE  FILTERS 

The  basic  criterion  of  performance  used  is  the  S/C  (signal-to-clutter) 
ratio.  Specifically  we  have  set  the  input  or  unfiltered  S/C  ratio  to  one; 
hence  our  S/C  ratio  is  actually  the  ratio  of  the  output  to  the  input  S/C  ratios. 
Some  authors  term  this  the  MTI  gain  or  improvement  factor.  The  loop  noise  is 
not  included  in  these  calculations;  however,  it  does  not  significantly  affect 
them. 

Another  important  criterion  is  the  rate  of  convergence  of  the  adaptive 
loops.  To  characterize  this  factor  we  h~ve  investigated  the  mean  S/C  ratio 
as  a  function  of  the  number  of  independent  samples,  subject  to  the  constraint 
that  the  time  constant  of  the  loops  is  chosen  to  give  a  constant  loop  noise 
factor  (n  =  0.1).  The  expected  actual  output  under  clutter  conditions  con¬ 
sists  of  two  factors — one  due  to  clutter  and  the  other  due  to  random  varia¬ 
tions  in  the  filter  weights  (loop  noise).  The  noise  factor,  n,  is  the  ratio 
of  the  output  due  to  loop  noise  to  the  output  due  to  the  clutter  when  the 

system  has  reached  steady  state.  (For  mere  detail  see  Appendices  A  and  C.) 

* 

The  initial  value  of  the  weights  has  been  chosen  to  be  GS  .  This  is  indepen¬ 
dent  of  the  clutter  spectrum,  optimum  when  only  receiver  noise  is  present, 
and  is  considered  a  reasonable  design  when  the  clutter  is  unknown. 

There  are  many  parameters  which  affect  the  performance  of  an  adaptive 
system.  We  have  chosen  to  keep  some  fixed  (e.g.,  the  antenna  pattern)  and 
have  normalized  others  (see  below).  To  depict  the  effect  of  all  of  these 
parameters,  we  have  chosen  a  baseline  set  of  values  and  then  (for  the  most 
part)  investigated  the  effect  of  varying  just  one  of  these  parameters  at  one 


time. 
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The  radar  backscatter  from  rain  has  been  modeled  as  a  gaussian-shaped 

spectrum,  its  mean,  variance,  and  power  relative  to  the  ground  clutter  being 

parameters.  In  the  case  of  rain,  the  initial  weights  were  not  only  chosen 
* 

to  be  GS  ,  but  in  some  cases  set  equal  to  the  mean  steady-state  value  obtained 
under  the  corresponding  no-rain  condition.  This  latter  initial  value  has 
been  termed  "rain  onset."  Interestingly  (for  the  cases  considered),  only  a 
slight  improvement  in  the  transient  response  takes  place  by  using  the  latter 
initial  values. 

The  maximum  MT1  gain  shown  in  these  curves  is  roughly  60  dB.  This  gain 
is  limited  by  the  sidelobe  clutter  spectrum.  A  sidelobe  level  of  29.6  dB, 
with  a  Dolph-Tschebychef f  antenna  pattern,  was  assumed  in  the  model.  The 
corresponding  two-way  sidelobe  gain  is  roughly  60  dB  below  the  maximum  main 
beam  antenna  gain.  This  sidelobe  level  is  representative  of  existing  radar 
systems.  Also,  roughly  50  dB  of  MTI  gain  is  adequate  in  most  AMTI  radars. 

In  practice,  the  performance  of  an  AMTI  radar  would  typically  be  limited  to 
50  to  60  dB  MTI  gain  by  receiver  noise.  By  assuming  a  lower  sidelobe  level, 
better  MTI  gain  could  be  attributed  to  the  adaptive  filters.  However,  this 
would  not  be  representative  of  typical  radars. 

3.1  PARAMETERS 

V  ■  mean  rain  velocity  (relative  to  ground) 

K 

VT  ■  target  velocity  (relative  to  ground) 

Vp  «*  platform  velocity  (relative  to  ground) 

A  »  wavelength 

f  *  pulse  repetition  frequency. 
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3.1.1  Antenna  Pattern  (Remains  fixed  throughout  the  plots 

depicted  herein) 

A  Dolph-Tschebychef f  design  with  elements  spaced  at  X/2.  Identical 
receive  and  transmit  patterns  and  zero  gain  in  the  reverse  half-circle. 
NEL  =  _30  =  number  of  elements 

SDLB  =  29.6  dfi  =  ratio  of  main  beam  peak  to  any  sidelobe  peak 
BEAM (DEG)  =  4. 1  =  beamwidth' (in  degrees)  between  3  dB  points. 

3.1.2  Ground  Clutter  (Assumes  independent  gaussian  stationary 

spatially  homogeneous  scatterers — 
clutter  spectrum  being  caused  by  the 
moving  radar  platform) 

NSC  =  100  =  number  of  equally  spaced  scatterers  (over  180°)  assumed 
in  computing  the  covariance  M 

NP  =  2,5,10,20  «  number  of  pulses  in  the  filter 

ALFA  =  a  =  1,_2,10  =  normalized  platform  velocity  =  2VpMfr 

SCAN  =  =  0(11.25)90° ,45°  =  scan  angle  =  angle  between  platform 

velocity  and  beam  center  (see  Figure  3) 

3.1.3  Rain  Clutter  (A  gaussian-shaped  spectrum) 

RNM  =  y  =  0.1  =  normalized  mean  =  2V  /Xf 

K  IT 

RNS  =  a  =  0.01,0.05,0.1  =  normalized  standard  deviation 

RNP  =  p  =  1_,10  =  ratio  of  total  rain  power  to  total  clutter  power. 

The  total  input  power  (rain  +  ground)  is  normalized 
to  one. 

3.1.4  Loop  (Applebaum  Array — see  Figures  1  and  2) 

O  C. 

GAIN  =  G  =  10,10  ,10  =  loop  gain 

ETA  =  n  =  0.1  =  noise  factor  =  ratio  of  the  output  due  to  random 

variations  in  the  weights  to  the  output  due  to  the 
clutter  in  the  steady  state.  This  and  the  gain 
determine  the  time  constant  of  the  loops. 
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3.1.5  Signal  (Normalized  to  unit  power,  point  source) 

BETA  =  3  =  0(0. 1)0. 9, 0,4  =  normalized  target  velocity  (relative  to 

ground) 

=  2V_/Xf  modulo  1 
T  r 

GAMMA  =  y  =  (acos^+3)  modulo  1  =  design  normalized  doppler  frequency 

for  the  filter. 

3.2  TYPES  OF  PLOTS 

3.2.1  Transient  Response 

The  S/C  ratio  (or  MTI  gain)  is  depicted  vs.  the  number  of  independent 

*  * 

samples.  The  initial  value  of  the  weights  is  chosen  to  be  GS  (where  S 
is  the  steering  signal  which  is  independent  of  the  clutter),  which  is 
considered  a  reasonable  design  when  the  clutter  is  unknown.  The  response 
at  this  initial  value  is  shown  by  a  small  dash.  Both  the  steady-state  and 
optimum  performance  are  depicted  as  straight  lines.  Due  to  the  high  base¬ 
line  gain  of  10^  on  most  plots,  the  steady  state  and  optimum  merge  into  a 
single  line.  Since  the  noise  factor,  r|,  is  kept  fixed,  these  plots  show  the 
convergence  rate  of  the  loops  at  the  same  relative  loop  noise  level.  Please 
note  the  logarithmic  "number  of  samples"  scale. 

3.2.2  Transient  Response  (Rain  Onset) 

Same  as  "Transient  Response"  except  that  the  initial  value  of  the  weights 
is  chosen  to  be  the  mean  steady-state  value  for  the  corresponding  no-rain 
condition.  The  small  dash  indicates  the  response  at  this  initial  value. 
Another  small  dash  labeled  "NO  RAIN  SS"  is  the  steady-state  response  under 
the  corresponding  no-rain  condition. 
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3.2.3  Clutter  Spectrum 

The  "folded"  (due  to  the  pulse  repetition  frequency,  f  )  spectrum  of 
the  clutter  is  depicted.  When  rain  is  present,  it  includes  the  contribution 
to  the  spectrum  by  the  rain.  It  is  normalized  to  unity  (0  dB)  total  power. 

3.2.4  Spectra-Filter ,  Clutter 

Depicts  two  folded  spectra  on  the  same  grid.  One  is  the  spectral  re¬ 
sponse  of  the  filter  formed  by  the  weights  (easily  distinguishable  since  it 
has  one  less  peak  than  the  number  of  pulses).  It  is  normalized  to  0  dB  at 
the  desired  design  frequency  of  y.  The  other  is  the  clutter  spectrum  as  spe¬ 
cified  above. 

3.2.5  Optimum  Performance 

Depicts  the  optimum  S/C  ratio  (MTI  gain)  vs.  the  normalized  target 
velocity  (B) .  A  number  of  curves  are  shown  at  various  scan  angles,  ip. 

3.3  DISCUSSION  OF  PERFORMANCE 

Very  good  steady-state  performance  (S/C  ratio  in  the  50  to  60  dB  range) 

is  obtainable  for  all  but  extreme  parameter  values  (just  what  these  extreme 

values  are  is  discernible  from  the  detailed  discussion  below). 

Reasonable  convergence  rate  of  the  loops  is  also  obtainable  for  most 

3  4 

parameter  values,  i.e.,  in  about  10  to  10  samples  the  mean  response  is  with¬ 
in  a  few  decibels  of  the  steady  state  which  is  negligibly  close  to  the  optimum. 
This  is  for  a  noise  factor  of  0.1;  hence  the  additional  degradation  of  per¬ 
formance  due  to  random  fluctuations  in  the  weights  at  this  point  is  negligible. 

Since  the  mean  weights  during  adaptation  are  not  equal  to  the  steady- 
state  weights,  the  loops  form  a  biased  estimate  of  the  weights.  This  causes 
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rather  slow  convergence  In  some  cases.  The  staircase-like  convergence  of 
the  mean  response  (even  on  a  logarithmic  scale)  indicates  that  there  are  long 
intervals  of  sampling  when  no  significant  improvement  in  the  mean  response 
takes  place  (whether  the  variance  of  the  weights  at  this  time  is  being  sig¬ 
nificantly  reduced  has  not  yet  been  investigated).  For  example,  for  the 
baseline  case  (Figure  4)  from  100  to  1000  samples,  there  is  only  about  a 
2  dB  improvement. 

For  some  extreme  cases  (very  heavy  rain)  (Figure  44),  the  mean  response’ 
even  diverges  (after  almost  achieving  the  steady-state  value)  to  a  value 
even  less  than  the  initial  value  (after  about  30,000  samples)  before  again 
converging  to  the  steady-state.  This  occurs  when  the  average  magnitude  of 

the  initial  weights  differs  widely  from  the  magnitude  of  the  steady-state 

-1  * 

weights.  The  steady-state  weights,  M  S  ,  tend  to  be  large  when  the  covari¬ 
ance  matrix  contains  one  or  more  very  small  eigenvalues.  It  might  be  prefer¬ 
able  to  have  an  unbiased  estimator,  i.e.,  one  in  which  the  mean  magnitudes 
of  the  weights  would  always  be  comparable  to  the  steady-state  weights.  In 
this  case  the  effect  of  additional  samples  would  be  to  adjust  the  relative 
magnitudes  and  phases  of  the  weights,  which  would  improve  the  response. 

Such  estimators  can  be  found — the  problem  is  that  they  markedly  increase 
the  complexity  of  the  calculations. 

It  is  felt  that  better  estimators  than  the  investigated  loops  are  pos¬ 
sible;  however,  probably  at  an  increased  complexity.  The  investigated  loops 
do  however  give  a  sufficiently  good  performance  that  is  probably  satisfactory 
for  many  applications.  The  simplicity  of  the  loops  then  makes  them  a  good 
overall  design  for  some  systems. 


3.3.1  Base  Line 


The  base  line  has  those  parameter  values  which  are  underlined  in  Sec¬ 
tion  3.  Figure  4  shows  the  mean  transient  response  with  its  characteristic 

staircase-like  (even  with  a  logarithmic  abscissa)  approach  to  the  steady- 

4 

state  value.  About  10  samples  are  seen  to  be  required  to  complete  adapta¬ 
tion.  The  duration  of  one  sample  is  the  radar  (compressed)  pulse  length, 
typically  around  1  microsecond. 

Figures  5  through  10  show  the  filter  response  (of  the  mean  value  of  the 
weights)  as  adaptation  takes  place.  The  clutter  spectrum  is  superimposed  on 
all  of  these  figures.  Initially  (samples  =  0)  no  specific  attempt  at  sup¬ 
pressing  the  broad  peak  of  the  clutter  spectrum  is  evident.  As  the  number  of 
samples  increases,  two. main  effects  are  discernible.  The  lobe  peaks  of  the 
filter  in  the  region  of  the  clutter  spectrum  maximum  are  decreasing  and  more 
filter  lobes  are  being  squeezed  into  the  extent  of  the  clutter  lobe  (from  2 
to  4  lobes).  The  main  filter  lobe  which  includes  the  selected  target  doppler 
frequency  also  widens  a  bit.  This  is  believed  to  be  a  necessary  consequence 
of  narrowing  the  other  lobes.  A  very  slight  shift  in  the  peak  of  the  main 
filter  lobe  may  also  be  noted.  This  effect  is  more  pronounced  under  more 
extreme  conditions  (e.g.,  see  Figure  51). 
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3.3.2  Optimum  Performance 

Figure  11  investigates  the  variation  of  performance  with  target  doppler 
velocity  (B)  and  scan  angle  Op).  The  optimum  S/C  ratio  is  shown.  However, 
for  a  loop  gain  of  10^  this  is  practically  (within  0.1  dB)  the  same  as  the 
steady-state  performance. 

As  expected,  performance  drops  sharply  as  the  target  doppler  approaches 
0  (or  1)  since  in  this  limiting  case  the  target  is  spatially  stationary  (or 
appears  so  due  to  the  doppler  ambiguity).  It  is  interesting  to  note  that 
performance  in  this  case  is  worst  at  0  scan  angle.  This  occurs  because  the 
platform  motion  clut ter  spectrum  becomes  more  spread  out  with  increasing  scan 
angle  (see  Figures  12  through  16),  thus  making  discrimination  possible. 
However,  as  the  target  doppler  moves  away  from  0  (or  1),  this  same  increased 
clutter  spread  causes  decreased  performance  since  it  is  more  difficult  to 
place  a  broad  null  in  tne  filter  spectrum.  The  erratic  behavior  of  the  per¬ 
formance  with  scan  angle  in  the  target  doppler  midrange  is  believed  to  be 
due  to  the  relative  difficulty  of  achieving  specific  separations  between  the 
filter  main  peak  and  predominant  broad  null.  In  the  midrange  of  target 
doppler  the  adapted  performance  is  seen  to  be  relatively  insensitive  to  scan 
angle,  showing  that  the  adaptation  manages  to  compensate  for  the  disparate 
clutter  spectra  (Figures  12  to  16). 


Figure  11 
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3.3.3  Variation  with  Scan  Angle 

3. 3. 3.1  Clutter  Spectra.  Figures  12  through  16  depict  the  variation 
of  the  clutter  spectrum  with  scan  angle.  Two  main  effects  are  to  be  noted. 
The  primary  peak  broadens  with  increasing  scan  angle  and  the  location  of 
this  peak  moves  towards  decreasing  normalized  doppler  frequency,  making  two 
complete  circuits  for  the  90°  variation  in  the  scan  angle.  These  two  cir¬ 
cuits  are  specifically  due  to  a  being  2.  The  spike  at  0  (or  1)  is  due  to  a 
being  an  integer  which  is  also  the  cause  of  the  spike  and  the  main  lobe 
being  superimposed  at  0  scan  angle  (refer  to  Appendix  A,  Equation  7). 

3. 3. 3. 2  Transient  Response.  Figures  17,  18,  4,  and  19  show  the  tran¬ 
sient  response  variation  with  scan  angle.  Note  that  at  0°  scan  angle  the 
initial  filter  does  extremely  well  since  the  clutter  spectrum  has  only  a  very 
narrow  peak;  however,  considerable  samples  are  required  to  improve  upon  it 
since  the  loops  find  it  difficult  to  place  a  very  deep  null  at  this  specific 
point.  However,  even  at  a  slight  scan  angle  (11.25°),  the  initial  filter  is 
seen  to  be  rather  poor  and  rapid  improvement  ensues  under  adaptation.  At  90° 
the  initial  response  is  even  worse,  yet  the  adapted  response  is  only  slightly 
worse  that  at  0°  scan  angle.  With  increasing  scan  angle,  the  transient 
response  becomes  more  and  more  step-like. 
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3.3.4  Variation  with  Target  Velocity  (g) 

(Figures  20,  21,  4,  22,  and  23)  It  is  interesting  to  note  that, 

though  the  steady-state  performance  changes  significantly  with  target  doppler 

(as  discussed  above),  the  number  of  samples  required  to  achieve  this  steady- 

4 

state  response  remains  essentially  10  samples.  A  more  and  more  pronounced 
staircase-like  transient  response  is  evident  with  worsening  steady-state 
performance. 
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3.3.5  Variation  with  Number  of  Pulses 

(Figures  24,  25,  4,  and  26)  With  increasing  number  of  pulses,  the 
variety  of  filters  possible  increases  rapidly,  so  it  is  not  surprising  that 
the  steady-state  performance  improves.  However,  the  complexity  of  the  system 
also  increases  markedly  so  that  by  20  pulses  the  additional  gain  in  performance 
probably  docs  not  warrant  the  increased  complexity. 
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3.3.6  Variation  with  Platform  Velocity  (a) 

(Figures  27,  4,  and  28)  Small  a  is  seen  to  give  better  steady-state 

4 

performance.  However,  it  is  to  be  noted  that  after  10  samples,  essentially 
the  same  performance  results  for  a  =  1  and  a  =  2.  Since  a  =  small  ct 

requires  a  low  platform  velocity,  high  pulse  repetition  rate,  or  a  long 
wavelength.  These  parameters  are  usually  decided  on  the  basis  of  other 
'system  requirements  than  filter  performance,  e.g. ,  the  PRF  may  be  selected 
to  avoid  second-time-around  returns  and  range  ambiguities. 
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3.3.7  Variation  with  Loop  Gain 

(Figures  29,  30,  and  4)  Note  that  the  time  constant  of  the  loops  is 
changed  so  as  to  maintain  a  constant  noise  factor  (n)  of  0.1.  It  is  seen 
that  low  loop  gain  causes  a  steady-state  performance  considerably  less  than 
the  optimum  performance.  The  transient  response  is  identical  for  the  various 
gains  except  that  with  lower  gain  the  asymptote  is  reached  sooner.  At  low 
loop  gains,  the  weights  approach  values  ([i+(I/G)]  S  )  which  are  different 
than  the  optimum  (M  S  ).  This  bias  in  the  steady-state  solution  is  negligi¬ 
ble  when  1/G  is  small  compared  to  the  smallest  eigenvalue  of  M. 
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3.3.8  Rain  Base  Line 

(Figures  31  through  36)  A  gaussian-shaped  rain  spectrum  with  mean  of 
0.1,  standard  deviation  of  0.05,  and  with  total  power  equal  to  the  total 
clutter  power  is  added  to  the  baseline  clutter  spectrum.  Note  that  the  total 
power  is  then  normalized  to  one  so  that  the  S/C  ratio  depicted  is  still  the 
ratio  of  output  to  input  S/C  ratio  (MTI  gain).  In  practice,  there  would  be 
an  increase  in  total  clutter  power  so  that  the  input  S/C  ratio  would  decrease. 
Thus  to  maintain  the  output  S/C  ratio,  the  MTI  gain  would  actually  have  to 
increase.  Specifically,  with  a  relative  rain  power  of  one,  the  MTI  gain 
would  have  to  increase  3  dB  to  main;ain  the  output  S/C  ratio  at  the  pre-rain 
level. 

Some  improvement  (about  4  dB  maximum)  in  transient  response  results  from 
initializing  the  weights  at  the  pre-rain  steady-state  values  (Figure  32). 
However,  as  the  number  of  samples  increases,  this  improvement  diminishes  and 
disappears  before  the  steady  state  is  approached.  The  baseline  rain  is 
rather  severe  (compare  Figures  5  and  33);  hence  the  adaptation  is  governed 
mostly  by  the  rain — not  the  ground  clutter. 

Note  in  particular  that  even  for  this  relatively  severe  rain  the  adapted 
MTI  gain  is  essentially  the  same  as  for  the  no-rain  case.  As  shown  below, 
for  lighter  rains  the  MTI  gain  may  actually  improve. 

Figures  33  through  36  show  the  filter  spectra  during  adaptation;  the 
clutter  (rain  and  ground)  spectrum  is  superimposed  to  show  how  the  filter 
nulls  out  the  broad  clutter  spectrum  peak. 
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3.3.9  Rain  at  0°  Scan  Angle 

(Figures  37  through  40)  The  broadness  of  the  rain  spectrum  (compare 
Figures  12  and  38)  can  clearly  be  seen  in  this  case.  Of  particular  interest 
is  the  tremendous  drop  in  performance  when  the  initial  filter  is  used  (recall 
that  in  the  no-rain  case  at  0s  scan  angle,  the  initial  filter  was  very  good, 
as  shown  in  Figure  17).  There  is  a  slight  (2  dB)  drop  in  the  MT1  gain  at 
steady  state. 
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3.3.10  Increased  Rain  Power  .  ? 

(Figures  31,  37,  41,  and  42)  Increasing  the  rain  power  relative  to 
the  clutter  power  has  the  curious  effect  of  increasing  the  adapted  MTI  gain. 
Note  however  that  the  input  S/C  ratio  would  now  drop  11  dB;  hence  the  output 
S/C  ratio  would  actually  be  worse.  Figures  31  and  41  show  the  comparison  at 
45°  scan  angle;  Figures  37  and  42  illustrate  the  effect  of  increasing  rain 
clutter  power  at  0°  scan  angle. 
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3.3.11  Variation  with  Rain  Spectral  Width 

(Figures  31  and  43  through  51)  When  the  rain  spectral  width  is  reduced 
to  0.01,  the  MTI  gain  increases  by  about  2  dB  over  the  no-rain  case,  which  is 
almost  enough  to  compensate  the  3  dB  drop  in  input  S/C  ratio.  Thus  the  effect 
of  this  narrow  spectrum  rain  is  almost  compensated  for  by  adaptation  and  the 
actual  output  S/C  ratio  drops  only  1  dB. 

Under  extreme  rain  spectral  width  (Rh’S  =  0.1),  a  very  curious  phenomenon 
occurs  in  the  mean  transient  response.  Adaptation  proceeds  reasonably  well 
until  4000  samples  (Figure  44),  but  thereafter  it  diverges  most  significantly 
before  converging  again.  Setting  the  initial  weights  to  their  no-rain  steady- 
state  mean  values  (Figure  45),  though  having  a  slightly  better  initial  per¬ 
formance  and  achieving  the  final  steady-state  value  a  little  sooner,  also  has 

a  more  serious  interim  drop  in  performance.  Of  particular  note  is  that  after 

* 

30,000  (after  initial  weights  of  GS  ),  the  performance  is  worse  than  initially 
Figures  46  through  51  show  the  filter  response  during  adaptation  for 
this  extreme  rain  case.  The  dip  in  performance  is  seen  to  be  due  to  a  large 
filter  lobe  at  other  than  the  design  frequency,  y This  is  believed  to  be 
caused  by  the  relative  changes  caused  in  the  transformed  weights  (see  Appendix 
A)  as  each  transformed  weight  in  succession  (determined  by  its  corresponding 
eigenvalue)  approaches  its  asymptotic  value.  The  latter  is  also  believed  to 
be  the  cause  of  the  staircase-like  transient  response.  We  have  not  yet 
found  a  simple  way  of  compensating  for  this  effect. 
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4.  SIMPLIFICATION  OF  ADAPTIVE  FILTERS 

The  complexity  of  adaptive  filters  which  are  implemented  digitally 

depends  strongly  on  the  number  of  bits  required  at  various  points  in  the 

system.  Both  the  data  storage  requirements  and  the  amount  of.  computation 

depend  on  number  of  bits.  Under  earlier  contracts,  TSC  has  investigated 

the  effects  of  quantization  noise  at  various  points  in  the  adaptive  control 

loops.  These  studies  were  performed  for  adaptive  array  antennas,  but  the 

results  apply  directly  to  the  closely  analogous  case  of  adaptive  filters. 

It  was  found  by  introducing  the  effects  of  quantization  in  an  adaptive  array 

simulation  that  a  single  bit  (one  bit  in-phase,  one  bit  quadrature)  is  suffi- 
* 

cient  at  the  v^  input  to  the  cross-correlator  (see  Figure  2).  An  analysis 
of  this  simplification  was  performed  under  this  contract  which  verifies  the 
simulation  results.  This  analysis  appeared  in  the  July  1972  issue  of  IEEE 
Trans.  AES,  and  was  included  in  the  interim  report^  on  this  contract. 

It  is  shown  that  the  estimated  covariance  matrix  differs  only  by  a  scale 

* 

factor  when  this  v^  input  is  represented  by  a  single  bit  per  quadrature 

component.  This  result  permits  a  major  simplification  of  adaptive  filters. 

Further  study  of  this  problem  has  shown  that  a  single-bit  representation 
* 

of  the  v^  input  to  the  cross-correlator  can  be  used,  i.e.,  one-bit  represents 
tion  of  just  one  quadrature  component.  The  analysis  supporting  this  conclu¬ 
sion  is  contained  in  Appendix  B. 


75 


5.  SUGGESTED  AREAS  OF  RESEARCH  AND  CONCLUSIONS 

As  discussed  in  Section  3  of  this  report,  the  Applebaum  type  of  adap¬ 
tive  array  gives  very  good  performance  for  its  relatively  simple  implementa¬ 
tion.  One  way  to  obtain  considerable  improved  performance  would  be  to 
estimate  the  covariance  matrix  directly,  to  invert  it,  and  to  calculate  the 
weights  by  W  =  M  S  .  This  procedure  is  not  as  complex  as  it  first  appears 
since  the  estimation  and  inversion  only  have  to  be  done  once  for  the  whole 
bank  of  filters  (each  filter  requires  its  own  set  of  loops),  the  different 
filters  being  obtained  by  a  matrix  multiplication.  The  inversion  of  the 
matrix  would  be  the  most  complicated  portion  of  this  procedure.  However, 

the  covariance  matrix  has  a  peculiar  nature  (i.e.,  M  is  a  function  only  of 

mn 

m-n — see  Appendix  A),  which  allows  it  to  be  inverted  in  a  particularly 
simple  manner  (see  [6]  and  [7]).  The  promise  of  a  considerable  improvement 
in  performance  at  a  modest  increase  in  complexity  makes  this  procedure  worthy 
of  investigation..  Such  an  investigation  may  also  shed  some  light  on  what 
minor  (in  complexity)  modifications  of  adaptive  loops  would  yield  the  most 
improvement  in  performance. 

The  sometimes  aberrant  behavior  (i.e.,  temporary  divergence  of  the  S/C 
ratio  during  adaptation)  of  the  loops  considered  does  raise  the  question  of 
whether  variations  of  adaptive  loops  (e.g.  ,  the  Widrow  array)  also  exhibit 
this  behavior. 

The  effect  of  the  loop  noise  on  the  S/C  ratio  has  been  determined  for 
the  steady-state  condition;  however,  its  transient  variation  during  adapta¬ 
tion  remains  to  be  investigated. 


76 


An  important  area  for  further  research  is  AMTI  radars  with  adaptivity 
in  both  space  and  time,  i.e.,  adaptive  control  of  both  the  antenna  pattern 
and  dopplor  filter  response.  One  important  optimization  study  in  space- 
time  adaptive  systems  is  related  to  the  choice  of  degrees  of  freedom.  For 
example,  a  system  with  12  degrees  of  freedom  could  use  a  single  antenna 
output  with  12-pulse  adaptive  filtering,  three  separate  antenna  outputs  on 
four  consecutive  pulses,  etc.  Antenna  outputs  could  be  obtained  from  sub¬ 
arrays  or  separate  beams,  e.g.,  sum  and  difference  beams  in  a  reflector 
antenna  with  multiple  feeds.  Methods  of  speeding  convergence  in  these 
systems  and  of  simplifying  the  implementation  are  important  areas  for  study. 

The  results  reported  here  show  that  adaptive  filters  can  provide  im¬ 
portant  improvements  in  AMTI  radars.  One  advantage  of  adaptive  systems  is 
the  ability  to  sense  the  presence,  of  rain  clutter  and  re-optimize  the  filter 

response.  The  convergence  rates  of  the  adaptive  filters  are  adequate  for 

4 

some  applications — the  typical  convergence  times  of  1000  to  10  samples 
correspond  to  1  to  10  milliseconds  in  a  radar  with  a  1  microsecond  compressed 
pulse  length.  The  change  in  clutter  due  to  antenna  scan  is  usually  small  in 
this  time  interval.  In  other  cases,  more  rapid  convergence  may  be  required, 
e.g.,  to  follow  changes  in  clutter  spectrum  with  range. 
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APPENDIX  A 


COMPUTATIONAL  EQUATIONS 

1.0  INTRODUCTION 

The  equations  used  in  computing  the  clutter  spectrum  (ground  and  rain), 
.covariance  matrix,  filter  weights  (optimum,  steady-state,  and  mean  transient), 
and  the  resulting  S/C  ratio  are  outlined  in  this  Appendix.  We  consider  the 
azimuthal  antenna  pattern  with  sidelobes  and  homogeneous  spatially  stationary 
clutter  Loth  in  azimuth  and  range.  We  neglect  considerations  of  the  varia¬ 
tion  of  antenna  pattern  due  to  radome  and  near  field  scatterers,  clutter 
inhomogeneities  due  to  terrain  type  and  incidence  angle,  elevation  antenna 
pattern,  and  depression  angle.  The  coordinate  system  ir  shown  in  Figure  3. 


2.0  NORMALIZED  PARAMETERS 

It. is  convenient  to  measure  the  platform  velocity  Vp,  and  the  target 
radial  velocity  (relative  to  ground),  V^,  in  terms  of  the  wavelength,  X,  and 
the  pulse  repetition  frequency,  f^. 

a  -2  Vp/Xfr 

B  -  2  VT/Xfr  (mod  1) 
y  *  a  cos(i{i)  +  g  (mod  1) 

Y  is  the  target  doppler  frequency  modulo  the  pulse  repetition  frequency. 
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3.0  EQUATIONS 

3.1  ANTHCNA  GAIN 

A  Dolph-Tschebyc.heff  antenna  pattern  (see  Ref.  3)  is  used.  The 
two-ray  antenna  power  pattern  is  given  by  (identical  receive  and  transmit 
patterns  are  assumed): 

-  G4«)  -  |  tk  [  Zo  J4  -f  <*  if 

=  0  othei*vise  (1) 

where 

T,.  =  order  Tschebycheff  polynomial 
li 

M  =  KEL-1  =29 

KEL  =  30  **  number  of  elements 

d  «*  A/2  =  spacing  between  elements 

-  1.01  ®  parameter  determining  sidelobe  level  and  beamwidth 

All  cases  considered  herein  use  the  same  antenna  pattern  with  parameter 
values  as  shown  above,  which  yield  a  main  beam  3  dB  width  of  4.1  degrees  and 
a  sidelobe  level  (ratio  of  r.ainbeam  peak  to  any  sidelobe  peak)  of  29.6  dB. 
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Vo  assume  spatially  stationary  homogeneous  clutter.  The  dopplor 
ccuer.cy  due  to  a  stationary  scattorer  at  angle  6  relative  to  the  platform 
iccity  (Fig.  3)  i s' 


f ,  =  2  V  cose/X  .  (2) 

d  p 

However,  the  dcppler  frequency  "folds  over"  at  the  pulse  repetition 
cqucncy,  f  ,  hence 

2  V  cosO 

fj/f  -  — -  =  a  cos0  (modulo  1).  (3) 

r 

The  clutter  spectral  density  is  thus 

S(f )  df  =  ^  G4  (0-i(O  |^||  df  .  (4) 

8  3  mod(a  cos0,l)  ®  f 


t  n  be  an  integer,  Arccos  be  the  principal  [0,ir]  inverse  function.  Then 


S(f)  - 


J2  g4w> 


d<(>  | 
df 1 


4-  =  Arccos 


df  +  2  ^ 

s  ~  ^  ♦  “  Arccos  j  +  if> 


(5) 


-sini^  <  <  1 

—  a 


sin4  <  ^  <  1 

—  a  — 


4 


Let 


Then 


sinO  =  cos^ 


T2  sini^ 


(for  the  1st 


summation) 


-  ?2  cos  +  T2  rin'J’ 


(for  the  2nc* 


summation) 


(6) 


Noting  that  =  1/oT^  for  either  summation  and  recalling  that 

G  (v)  =  g(s:in£)  v/e  obtain 


S(f)  =  z  g(T  cosiJj  -  T  sini|)/uT^  +  g(“T  cos’il  ~  T2  sinilO/aTj 

ntN^  neN^ 

vhcre 


(7) 


g(sin^)  is  defined  in  Eq.  1. 
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3.3  cKor:n  clutter  covaaIAkce 

Let  N  be  the  number  of  pulses  to  be  processed  coherently.  The  n^1 
pulse  voltage  return  from  a  scatterer  at  angle  0  Is  given  by  (normalized 
to  the  first  pulse): 


a  =  e  12™<n-,)  cos0g2(6-;-) 
n 


(8) 


The  contribution  of  this  clutter  element  to  the  clutter  covariance  is  given  by 


*  T 

MG=A  a 


(9) 


where  A  is  column  vector  of  the  a  . 

n 


Thus  f 


•G  (n,n)  = J" ( 


-i2u(m-l)acos9  „ 

;  GZ(6-i^) 

i2tf(n-l)acos9  _ 


x  e 


G  (G-ifO  d9 


■/« 


-i2ir(m-n)acos9 


G  (9-iJi)  e 


d9 


(10) 


Since  G  ($)  a  G  (— <f>) ,  G($)  =  0  for  |<j>|  >  ir/2 


Mc(n,n) 


r/2  . 

- /c*<*)  [ 


-i2tr(m-n)acos((})-0  -i27r(in-n)acos((f>+\j,)J 


+  e 


d$  (11) 
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The  preceding  integral  is  approximated  by  a  sum  of  NSC/2  equally  spaced 
angles  ever  the  region  (0,n/2).  This  matrix  is  scaled  so  that  the  diagonal 
terms  (which  are  all  equal)  have  the  value  of  which  corresponds  to  unity 
input  clutter  power.  Note  that  the  covariance  matrix  is  really  only  a 
function  of  n-n. 

3.  A  P.AIN  C!  UlTr.R  STTCTid’H 

Following  Nnthanson  [Ref.  4],  we  assume  that  the  rain  spectrum  is 
Gaussian.  It  is  again  convenient  to  introduce  normalized  parameters  for 
the  mean  rain  velocity  (relative  to  ground),  V  ,  and  the  standard  deviation 
of  the  spectrum  (mainly  due  to  wind  shear  effects),  oR.  Thus  let 

V  RNM  =  2  VD/Xf 
R  r 

a  =  RNS  =  2  oD/Xf  . 

R  r 

Due  to  platform  motion,  the  normalized  mean  relative  to  the  radar  is 
p  +  acosijj.  The  rain  clutter  spectrum  (normalized  to  unit  power)  is  hence 

-  (f-u-cicosi|')^/2a^ 

S -(f)  -  — “  e  (12) 

/2T  a 

The  above  analysis  ignores  the  "fold  over"  of  the  doppler  frequency  due 
to  the  pulse  repetition  frequency,  however,  it  is  an  excellent  approxima¬ 
tion  as  long  as  o  is  say  less  than  0.1. 
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3.5  FAT*;  CLUTTER  COVARIANCE  ' 

The  rain  clutter  correlation  function  is  the  inverse  Fourier  trans¬ 
form  of  its  spectrum,  i.e, 

R(t)  -^_1SR<f).  =  /  SR(f)e32llfTdf  .  (13) 


From  Reference  5 


her.ee. 


j2ir(y+acosi{j)f  t  —2 (Traf  t)^ 
R(t)  -  e  r  e  r 


(14) 


and  the  clutter  covariance  matrix  is 

2 

j2r(vi-Ktcos^)  (m-n)  -2(Tto(ra-n)) 

^(m.n)  =  e  e  (15) 

3.6  COMBINING  RAIN  AIT)  GROUND  CLUTTER 

Let  p  =  RNP  be  the  ratio  of  the  total  rain  power  to  the  total  ground 
power.  The  rain  and  ground  covariances  and  spectra  have  been  normalized 
to  unit  power.  Hence,  since  the  rain  and  ground  clutter  are  assumed 
statistically  independent,  the  total  covariance  (also  normalized  to  unit 
power)  is  given  by 
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M(m,n)  =  MG (13,11)  +  M^Gn.n)  (16) 

Mote  again  that  all  the  covariances  are  really  only  functions  of  m-n. 
The  total  spectrum  is  given  by 

S(f)  ‘  W  SG<f)  +  SK<f>  •  <17> 

3.7  OPTIMUM  FILTER 

As  has  been  shown  in  Section  2  of  this  report,  the  weights  for  an 
optimum  filter  are  given  by  (or  proportional  to): 


Since  the  eigenvalues  and  eigenvectors  of  M  will  be  required  for  the 


analysis  of  the  transient  behavior  we  obtained  the  optimum  weights  by  a 


different  computational  procedure.  We  compute  Q  and  A  such  that  A  is 
*T 

diagonal,  Q  Q  =  I,  and 


Note  that  since  A  is  diagonal,  the  computation  of  A~*  is  trivial. 


9 


(24) 


10 


11 


G(At)Tr 

We  define  n  =  C  -  1  =  — — -  as  the  loop  noise  factor.  Then  the  mean 

transient  weights  are 


w  (i)  =  [  w 
n  1  no 


X  +  1/G 
n 


-a  +  i/G) 
n  i 

e  r 


X  +  1/G 
n 


where  i  =  the  number  of  independent  samples,  spaced  At  apart 

G  =  control  loop  gain  (determined  by  allowable  steady-state 
degradation  from  optimum  performance) 

=  trace  of  covariance  matrix  (total  clutter  power  in  all  N 
pulses) 

n  =  loop  noise  factor  (was  shown  to  be  0.1) 

3.12  mean  tra::-ste::t  response 

We  have  calculated  the  mean  transient  response  using 

*  * 

w_  S  S-  v 

«/C  '  -5 - — 

wT  A  w 

where  w  are  the  mean  transient  transformed  weights.  This  neglects  the 
effect  of  loop  noise  on  the  S/C  ratio,  however,  loop  noise  is  by  design 
(choice  of  noise  factor  .)  quite  small  so  the  above  approximation  is  close. 
As  shown  in  Appendix  C,  the  2nd  order  transient  moment  of  the  weights 
is  required  to  determine  this  loop  noise,  which  can  be  determined  by  an 
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iterative  procedure.  For  the  steady-state  a  direct  solution  has  been  found 
(i.e.,  the  noise  factor)  but  no  simple  solution  for  the  transient  loop 
noise  has  yet  been  found. 


APPENDIX  B 


GENERAL  EFFECT  OF  ENVELOPE  NORMALIZATION  IN 
ADAPTIVE  ARRAY  CONTROL  LOOPS 

G.  W.  Lank 


introduction 

The  properties  of  an  adaptive  array  antenna,  including  transient 
response  rate  and  control  loop  noise,  depend  on  the  intensity  of  the  exter¬ 
nal  noise  field.  The  dependence  can  be  reduced  by  a  general  envelope 
normalization.  This  can  be  done  without  degrading  the  performance  of  the 
adaptive  array.  Special  cases  of  the  normalization  are  envelope  limiting, 
considered  in  [4],  and  one-bit  digitization  of  the  real  and  imaginary  parts 
of  the  signal  from  which  the  envelopes  are  formed,  considered  in  [5], 
Another  important  special  case  is  one-bit  digitization  of  the  imaginary 
part  of  the  signal  from  which  the  envelopes  are  formed  while  the  real  part 
is  set  to  zero. 

DISCUSSION 

Adaptive  array  antennas  have  been  discussed  in  [1],  {2],  and  [3]  by, 
respectively,  Widrow,  Applebaum,  and  Brennan,  et  al.  In  [4],  by  Brennan 
and  Reed,  it  is  shown  that  envelope  limiting  in  the  control  loops  reduces 
the  effects  of  varying  noise  intensity  without  degrading  array  performance. 
In  [5],  by  Lank  and  Brennan,  this  is  also  shown  to  be  true  when  one-bit 
digitization  of  the  re.nl  and  imaginary  parts  of  the  signal  from  which  the 
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envoi ope  is  obtained  is  performed  instead  of  envelope  limiting.  It  will  be 


shown  that  the.  results  of  [4]  and  [5]  are  special  cases  of  a  general  en¬ 
velope  normalisation  to  be  defined. 

Another  important  special  case  of  the  general  result  is  to  digitize  to 
one  bit  the  imaginary  part  of  the  signal  from  which  the  envelope  is  obtained 
while  setting  the  real  part  to  zero.  This  eliminates  one-half  the  multipli¬ 
cations  in  the  control  loop  (as  well  as  eliminating  the  storage  of  the  real 
part  of  the  signal  from  wh.  ,h  the  envelope  is  formed  in  some  adaptive  systems), 

The  loops  considered  are  those  analyzed  in  [3]  and  [4]. 

til  & 

Let  vjn  =  the  complex  video  input  to  the  m  array  element,  v  =  the 

* 

complex  conjugate  of  v^,  and  f (v  )  =  the  function  of  the  envelope  applied  to 
the  control  loops. 

As  a  direct  consequence  of  the  results  of  [3],  the  average  value  of  the 
transient  response  of  the  loops  as  well  as  the  steady-state  PUS  noise  are 
only  dependent  on  the  covariance  matrix 


A  *=  E[f(v  )v  ] 
mn  m  n 


This  is  true,  assuming  the  loop  time  constants  are  long  compared  to  the 

correlation  times  of  v  and  the  steady-state  RMS  noise  is  small. 

n 


v  °Rc 
m  m 


Then  assume 


f(vn>  "  f(Rnc  "  c(e 


Thus  (3)  nnsur.es  that  the  function  of  the  envelope  applied  to  the 

loops  is  only  a  function  of  <{>  .  Hence  it  is  a  general  normalization  in 

that  f (v  )  is  independent  of  R  .  ■' 

ei  m 

References  [4]  and  [5]  can  be  considered’ special  cases  of  (3).  In  [4], 


.  I  -W„\  -J*. 


f(vn)  »  E 


n  - 


while  in  [3], 


l(\?  “  «(e  J  i  = 


1-j  „  OS  <f>  <  -rr/2 

-1-j  ,  71/2  S  <f>m  <  IT 

-1+j  ,  ir  ^  ♦m  <  3ir/2 

1+j  ,  3rr/2  S  $  <  2n 


As  in  [4]  and  [5],  we  shall  calculate  the  covariance  matrix  in  (1)  for 
the  case  where  is  a  zero  mean  complex  stationary  Gaussian  process. 

The  joint  probability  density  of  R^,  R^,  is  given  as  in  [41  by 


R  R  f  R2+R2-2pR  R  cos  (6  -6  -!&) 

m  n  1  m  n  in  n  n  m 

p  (R  ,K  ,9  ,9  )  K  — r- 7 - tt~  exp  <-  - = - - 

r  ra’  non  .  2  4.,  2.  r  }  _  2.,  2. 

4«  o  (l-p  )  (  2a  (l-p  ) 


where 


°2  "  I  ^{hml2}  *  i  E{lv„|2}  l  E{vvj  "  -  H_ 

E{Vn}"° 
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the  above,  the  elements  of  the  covariance  matrix  A  arc 


A 


2  v  co  2tr  .  ^ 

f  dA  J  dR  f dR  f  d'-I>  o/e  3  n)R  e3  np(R  }>  4,  <f>  ) 

J  sn  J  m  j  nj  n°  \  /  11  m’  n*  m*  n 


(7) 


l,ct  f,  =  (J>  -  s'  in  (7).  Then,  using  (6)  and  the  periodicity  of  the 

cosine  function,  (7)  becomes 


inn 


mh(  1,»E(e  *°) 


Mu 

e  d<J» 


m 


(8) 


where 


in  co  Oy 

f  f  7 

J  -  /  dR  dR  f  <v2 
un  J  m  J  a  J 

0  0  0 


R  R2 
rn  n 


(  RZ4R-2pR  R  cos 3  ) 

- ^ - —  exp  j-  - - j - - -  +  J  C3  K0/ 

2ira  (1-p  )  (  2a  (1-p  )  J 


(9) 


In  [4],  p  is  obtained  (see  Eqs.  (8)  to  (13)  inclusive  of  [4]),  as 


; —  M 

p  -  — 

Mnn  \ 8  0 


(10) 


Thus,  substituting  (10)  into  (8),  one  obtains 


M 


rai 


ran 


o/32tT 


•  f  * 

•ti 


DI 


d«, 


ta 


(11) 


4s  in  [4]  and  [5],  the  dements  of  the  covariance  matrix  A  differ  from 
the  elements  of  M  only  by  a  common  factor  which  is  proportional  to  l/o.  Thus 
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APPENDIX  C 

KEAN  STEADY  STATE  OUTPUT  CLUTTER  -  APPLEBAUM  ARRAY 

In  an  adaptive  filter  or  array,  the  weights  are  a  stochastic  process. 
Expressions  for  the  steady  state  mean  and  covariance  of  the  weights  are 
derived.  From  these  an  explicit  relationship  for  the  mean  steady  state 
output  clutter  is  calculated.  An  iterative  solution  to  the  mean  transient 
output  clutter  is  inherent  in  the  method. 

1.0  System  Descrintion 
1.1  Schematic 


(50- 


A 


: *1 


w 


V  • 


. 4 


Solid  lines  denote  column  vectors 


Output 


Fig.  1  -  Applcbaum  Array  (Filter) 
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1.2  Eciuations 


iii+  [v*vT  +  i/g]w  »  s* 

G 


(i) 


V  (when  no  signal  is  present)  is  a  vector  of  zero  mean  complex 
Gaussian  stationary  random  variables  with  covariance 

EV*VT  =  M.  (2) 

where  E  denotes  the  expectation  or  average. 

Output  =  VJT*  V*VTW  (3) 

Note  that  EW^  V*VTW  4  EWT  EV*V^  EW  unless  the  tine  constant  of  the 
integrator,  x,  is  chosen  sufficiently  large. 

2.0  Difference  Equation 

Let  At  *»  time  between  samples  of  V,  then  (1)  becomes 

Vat  •  -f  [s*  -<vJvt  +  I/G)wt]+  <5> 

3.0  Steady-state  naan  weights 

Taking  the  expectation  of  (5),  noting  that  since  Vf  is  determi  cd  by  V 

t  u 

u<t-At  it  is  statistically  independent  of  V  , 
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“wit  ■  ^  [  s*  •  <«) 

In  the  steady  state  U  . ,  =  W  =  W  hence 

.  t-rAt  t 

(M  +  1/G)W  =  S* 

W  =  (M-rT/G)  1  S*  For . convergence:  (A/fl/G)  <1.  (7) 
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6 


(24) 
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Adaptive  filtering  is  a  technique  for  optimizing  the  doppler  filter  response 
in  an  Mil  (moving  target  indication)  radar.  In  aiiborne  MTI  radars,  the  clutter 
spectrum  is  continually  changing  and  varies  with  scan  angle,  antenna  pattern,  and 
angular  distribution  of  the  clutter  intensity.  Rain  backscatter,  if  present,  has 
a  doppler  spectrum  which  depends  on  local  wind  velocity  and  wind  shear.  An  adaptive 
filter  senses  each  of  these  effects  and  adaptively  controls  the  filter  weights  to 
maximize  the  signal-to-cl ut ter  ratio  in  the  filter  output.  Curves  are  presented 
in  this  report  which  illustrate  the  performance  of  adaptive  filters  in  AMT I  radar. 
The  steady-state  and  transient  response  of  adaptive  filters,  and  performance  in 
rain,  are  included.  Methods  of  simplifying  adaptive  filter  control  loops  are 
discussed. 
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